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I. INTRODUCTION 

In this paper, we derive simple new analytic expressions for a standard D dimensional 
one-loop two-point Feynman integral with two different masses > and m| > and an 
arbitrary external momentum p^, 



Our results are given terms of a single, manifestly symmetric Appell function^ F4, or alter- 
natively, the Appell function Fi. 

The integrals like that in (1) are ubiquitous in quantum field theory and renormalization 
group theory of phase transitions and critical phenomena. Actually, there exists for long time 
a standard general result for the integral (1) given by Berends, Davydychev, and Smirnov 
in Ref. 2 henceforth called BDS. This result was expressed essentially in a form of a linear 
combination of two Gauss hypergeometric functions of complicated arguments. 

Clearly, the Feynman integral (1) is a symmetric function of its parameters ml and ml. 
In the result of BDS, which we shall quote further, this property was implemented through 
a symmetrization 



of two complicated terms, none of which obeyed the symmetry mf <-> m|. 

Here we present the alternative general results for the integral (1). These are given 
in terms of a single Appell function, F4 or Fi, manifestly symmetric with respect to the 
interchange of mf and mg. The symmetry mf ^ m\ is implemented within the arguments 
of these functions. That is, our new expressions are given in terms of completely symmetric 
combinations like m\m\ or (mi ± 7712)^. They do not contain any symmetrizations like that 



We believe that such kind of results are interesting on their own right, and hope that 
their existence can be useful in the quantum field theory. Moreover, comparing our results 
with the expression of BDS we arrive at special identities between the Appell and Gauss 
hypergeometric functions that are apparently new mathematical relations. 

Moreover, there exists another reason for practical utility of our results. In fact, our 
motivation for present calculations was in the theory of the critical behavior of statistical 
mechanical systems at Lifshitz points^. Reviews and extensive lists of references on this 




(1) 



F[m^, 777.2) + ^{^21 ^1) 



(2) 



in (2). 
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topic can be found in Refs. 4-7. 

The field theoretical treatment of m-axial Lifshitz points in d-dimensional systems fre- 
quently requires calculations of momentum integrals 



in contexts of the ^-expansion (see e.g. Ref. 8), or large- A?" expansion^. In this setting, 
an m-dimensional subspace corresponding to < m < special anisotropy directions must 
be split out within the whole (i-dimensional space and complemented by the remaining 
(d— m)-dimensional subspace with no peculiarities in its physical properties. The momentum 
integral in (3) is over M*^ = MJ^ © ]R'^~'", and the inner integration over p in W^~'^ is easily 
recognized to be the same as in (1) with D = d — m and identifications of masses m\ and 
m| as (qr ± q^/'^Y. In the fully isotropic limit m = 0, we recover the usual situation of the 
Euclidean field theory in = 4 — £ dimensions with the standard upper critical dimension 
D* = 4. The latter follows from d*{m) = 4 + m/2, the upper critical dimension appropriate 
for the Lifshitz-point theory, as the number of the anisotropy axes m vanishes. 

In the generic case of non-isotropic m-axial Lifshitz points with < m < d, when we are 
interested, for example, in the one-loop energy-energy correlation function, a further integral 
over q has to be still performed in (3). To be feasible, its integrand must be as simple and 
symmetric as possible. We could not use the BDS result of the inner p integration to proceed. 
The reason was the lack of symmetry mf ^ m\ in individual terms F of (2). This implied 
the presence of complicated explicit dependencies on odd powers of the integration variable 
q in each of these terms. Of course, they cancel out in the symmetric sum (2), which is an 
even function of q (and qx). But the outer integration of any part of this sum appeared 
completely hopeless. 

In a hope to achieve the desired simple symmetry of the inner integral in (3), we have 
differently tried to apply the numerous standard transformation formulas of Gauss hyperge- 
omctric functions in the BDS result, but in vain. The application of the Feynman-parameter 
integration also lead us only to results equivalent to that of BDS. Hence, we were urged to 
attempt a completely different calculation in order to get the inner integral given separately 
by (1) in a simplest possible and manifestly symmetric form. 

Before turning to details of our calculation and its implications, we write down the BDS 




J{Px,qx) 
f d'^q 



1 1 



(3) 



^2n)d-m (p _ pj2y + {q- q,/2)^ {p + qj2y + {q + q,/2)^ 
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result in the following section. 



II. THE RESULT OF BERENDS, DAVYDYCHEV, AND SMIRNOV^ 

The Feynman integral of equation (1) obeys the scaling relation 

I{p^; ml ml) = 7(1; mj/pl, ml/pl) (4) 

where we define, as usual for the critical phenomena theory, D = 4 — e (note that this differs 
from the space-time dimension n — 4 — 2e ol BDS). Hence, without loss of generality we can 
consider the function (cf. (3)) 

I{ml ml) ^ /(I; ml ml) = J ^^^n _ yly + ^2 (p+ 1/2)2 + ^2 (5) 

where 1 denotes an arbitrary vector of unit length. The dependence of the original inte- 
gral (1) on external momentum can be easily reconstructed by scaling at any stage of 
calculations. 

The reference 2 contains the explicit result for this integral in terms of Gauss hyperge- 
ometric functions. Before quoting it, let us split from the function 7(m^,m|) an overall 
numeric factor via 

I{ml ml) = (47r)-f r(2 - f ) /(m?, mj) . (6) 

While r(2 — D/2) is singular as D ^ 4, the non-trivial function I^mlm^) is finite in this 
limit. In our conventions, the translation of the BDS result (see Eq. (A. 7) of Ref. 2) reads 

m?-\, 2 r.f. n D D {1 + ml - mj - VA)\ 



+ ^{l+ml-mi-^),F,[l,2-f,f, 



4mf / 

+ {ml ^ ml) . (7) 

For brevity, we denoted by (m^ ^ ml) the presence of the third term, repeating the function 
given explicitly in the second line, but with interchange of the masses mf and ml The 
symmetric combination 

A = {ml + ml + 1)^ - Amjml = [(mi + 7712)^ + l][(m2 - mi)^ + 1] (8) 

is a counterpart of the Kallen function with unit momentum in the terminology of Ref. 2. On 
the other hand, this is the discriminant of a quadratic equation that has to be solved in the 



course of calculations employing the Feynman parametrization. The apparent singularities 
at D — 2m (7) mutually cancel, and the finite result in this case is given explicitly in section 



The original BDS result was derived using the technique of Mellin-Barnes contour integral 
representations^". The result (6)- (7) or its immediate generalizations, expressed in terms of 
generalized hypergeometric functions of two variables, have been reproduced several times 
by different authors using different means^^"^^. All of them either contained symmetrizations 
like (2) or comprised some hidden symmetries in apparently non-symmetric expressions. 

As discussed in the Introduction, for our purposes we needed a completely different kind 
of result for the function /(m^,m2). Thus, wc had to calculate the integral (5) by using 
another procedure, not related to that of BDS or the Feynman parametrization. The way 
of doing it is described in the next section. 

III. THE ALTERNATIVE CALCULATION 

Let us return to the integral (5). Denoting the denominators of its integrand by A_ and 
A+, we write it as 



Using the partial fraction expansion we decompose our integral into two terms via 



Now, it is useful to introduce the arithmetic mean a — {A^ + A+)/2 and to express the 
denominators A_ and A^ in terms of their mean value a and deviation from it b, as A± — 
a±b. Hence, 



This representation brings into consideration the mass center — {ml + m\)/2 and the 
deviation = {m\ — m\)/2 ^ of masses m\ and m\ from m?. In terms of these values, 
we have 



V. 




with 



A_ =pHm2 + l/4-(p-l), 
=p2^mi+l/4+(p-l). 





a — -\- 
6 = (p • 1) + . 



with = + 1/4 and 
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The above simple algebraic manipulations provide us with a basis for further calculations 
in terms of certain combinations of ml and ml, more appropriate for the parametrization of 
final results than the original masses itself. These allow us to write 

I{mlml)^l{I_ + U) 



with 



{27t)^ + ml + ml ± mf, ± (p ■ 1) 
Here is a manifestly symmetric value related to the mass center m^ . On the other hand, 

the value mi, is given by the difference of m^ and m\. It is not symmetric with respect 
to the interchange of masses and thus still not "good". Besides of dependency on m^, the 
result must depend on m^, not m;,. Our aim is now to reach a manifest realization of the 
symmetry m\ <-> m\ at an early stage of the calculation, not on the level of final results. 
As a next step, we perform the angular integration in I± via 

^ ^ Jo p^ + ml Jq VL p^ + ml ±mb±p cos 6 

Here Kd is a usual geometric factor given by Kd — 2^~^7r~^/^/r(£)/2). The normalization 
factor of the angular integration fl — ^/t^V{D/2 — l/2)/T{D/2) is the value of the 9 integral 
without the denominator. The required result is 

pD-i I ^ ( I D p2 



r dp ^ ^ F (i 

Jq p^ + ml p^ + ml ± mfe V ' 2 ' 2 ' 



[p^ + ml ± mijY ) 

Hence we obtain, by definition of the Gauss hypergeometric function, a series represen- 
tation 

J , Kq ^ (l/2)„ /•» P°^'+^" 



X 



(l-^) +( 

\ p'^ -\- / \ 



2n-l / jYi^ N-2n-l 



+ 1 + 



pj. _|_ / J 

where (c)n = T[c + n)/T[c) is the Pochhammer symbol. Here, inside of the square brackets 
we get a simple realization of the symmetry m\ ^ m\ in the form of (2). This symmetry 
becomes manifest in the representation of the content of square brackets in terms of a Gauss 
hypergeometric function (Ref. 15, p. 34, Ref. 16, p. 461, Eq. 106) 

/ m,i, \ -2"-i / mh \ / I 1 ml \ , , 

+ l + ^T^ ^22Fi n+-,n+l;-; " 10 

\ p^ + mlJ \ p^ + m^/ \ 2 2 {p'' + ml)^J 



the argument of which depends exphcitly on ml — (m| — niD^/A. 

Using equation (10) we can simply perform the remaining integration over the radial 
distance p in (9). Its output is again (see (6)) given as I{ml,ml) = (47r)~^/^r(2 — 
D/2)i{m1,m2), while for the function I{mf,m2) we obtain a double series expansion 

The same result could be reached directly from equation (9) by using binomial expansions 
for each of two terms in the square brackets. 

The double series in (11) does not match directly any item of the Horn list (see Refs. 
15,17) of essentially distinct complete hypergcometric functions of two variables. This leaves 
a possibility that it can be reduced to some standard function through certain algebraic 
rearrangements . 

We see that one of the Pochhammer symbols in the nominator of (11) has a doubling of 
the summation index k. This is rather unusual in the theory of generalized hypergcometric 
functions. In order to reduce our double series to a more common form, one has to get rid 
of this duplication. This can be achieved by applying the exact resummation formula (Ref. 
15, p. 100) 

n 

J2 C(n, A;) = ^ ^ C7(n - fc. A;) 

k,n>Ci n>0 k=0 

along with the identity (Ref. 15, p. 23) 

1 (-l)^(-n)fc ^ , 

7 7TT - ^ ^ ' , 0<k<n. 

[n — ky. ni 

These manipulations lead us to 

Here we identify the inner finite sum with 2F'i{— n,e /2 + n] 1/2] —ml/ml) and convert it 
to a non-terminating hypergcometric function with the help of the linear transformation 
formula (see e.g. Ref. 18, p. 559, Eqs. 15.3.4-5) 

2Fi(a, 6; c; ^) = (1 - z)-'' 2F, (a, c - 6; c; . (13) 

Thus we get 

f. 2 2^ (i)n(|)n / 1 \" f £ , 1, 1 4mh 

H<ml)={-) ) ,F.(- + n,- + „;-;^) (14) 

n>0 ^2/" 



where we again encounter the famihar discriminant A from equation (8). This time it 
appears in a compact form as 

A = 4(m2 + mg). (15) 

The sum in equation (14) represents the Appell function F4 (Ref. 15, p. 53, Ref. 16, p. 
413, Eq. 5) whose original definition is^'^^'^^'^*^ 

^.(«,6;c,c';.,.)^i:^^|;f|flJS ^+/M<1. (16) 

k,n>0 \ ■> \ ■> 

Thus, the basic integral I{ml,ml) = (47r)~'°/^r(£/2) I{ml,ml) is given by 

ff 2 2n f^Y^f^ 13 11 4ml\ 

7(m„m,)=(-j ^4(-,-;-,-;^,^j. (17) 

For reader's convenience, we recall the notation = {ml — m\)/2 and note an alternative 
exphcit formula (cf. (8)) for the value A, which follows directly from (15): 

l^^l + 2{ml + ml) + {ml-mlf . (18) 

Let us check the convergence region of the Appell function F4 in (17). Following (16), we 
have to make sure that the combination 

1 I 2|mfe| _ l + |mi-m^| ^ ^^^^ 



a/A a/A 

To this end we rewrite the value A as 

A = {l+\ml- mWf + 2{ml + ml - \ml - m\\) 
= {1 + \ml — m\\Y + 2m.va.{m\,ml}. 

The last equality gives evidence that \/A > 1 + |m| — mf|, and thus, the inequality (19) 
holds for any m\, ml > 0. 

The equation (17) is one of our central results. The integral defined by (5) is expressed 
here in terms of a single Appell function F4. Its variables are given in terms of manifestly 
symmetric in mf and m| combinations A and mf. The result for 7(mf,m|) is finite and 
absolutely convergent in the range < D < 4 for arbitrary non-vanishing positive masses 
mf and ml- The limit D — > will be considered in section V along with further special 
cases. 
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IV. A FURTHER TRANSFORMATION OF THE APPELL FUNCTION 



There is a repetition of two parameters in the Appell function F4 appearing in the equation 
(17). This means that the result (17) can be alternatively expressed in terms of the Horn 
function (Ref. 15, p. 57). Moreover, it is reducible to the Appell function 15,19,20^ 

Fi{a;b,b;c;x,y)^ y — -— , max{|x|, |y|} < 1, (20) 

k,n>0 ^^^*^+" 

by means of Bailey's^^ reduction formula^^ 

-X -Y 



F4(a,b; c,b;—— 



x)(i-F)'(i-x)(i-r), 

= (l-X)"(l-F)"Fi(a;c-6, l + a-c;c;X,XF). (21) 

Denoting for a while the arguments 1/A and Ami/ A Appell function F4 in (17) as 

and B^, we find (cf. Refs. 22,23) 

2B^ ' 2A^ ^^^^ 

with the new discriminant 

5 = (A" + - If - AA'B^ = [{A + Bf - 1][{B - Af - 1], (23) 

astonishingly similar to the function A from equation (8). 

The last three equations yield totally unexpected simple expressions. In particular, we 

get 



v^ = 4 



2^2 



m 



2 



A 

Choosing the upper + sign in the solutions X and Y from (22) and assuming that mi > 
and m2 > are the principal square roots of m\ and m|, we obtain 

X^-- ^ and y = -(mi-m2)'. (24) 

[nil + "^2) 

An alternative choice of minus signs in (22), as in Ref. 22, gives X — — l/(mi — 777.2)^ 
and Y — —(mi + m2)^. This would lead to a much less convenient parametrization of the 
function I{m\, ml) owing to the singular behavior of the variable X by equal parameters 
mi and m2. 
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Thus, the apphcation of the reduction formula (21) to the function F4 from (17) leads to 
the Appell function Fi with surprisingly simple symmetric arguments. For I{'ml,'ml) this 
yields 

mi + m2\-'^ f e 1 e 3 -1 {nii - 



Reintroducing by scaling (4) the momentum dependence of the original Feynman integral 
(1) and taking into account the overall numeric factor (6) we obtain a very appealing in its 
structure result 

7(p,;m?,m^) = (47r)-^+^/2r(|) (26) 



n^^j ^^U'^'-2 + 2'2' 



-pI (mi - ms)^ 



(mi + ' (mi + m2)^ , 

A note is in order here. The last two equations represent the double series expansions 
(20) whose domain of absolute convergence is shrunk with respect to the original integral's 
validity range, mi > 0, m2 > 0, and arbitrary p^. Indeed, the definition (20) of the Appell 
function Fi requires mi + m2 > 1 for the first variable of the scaling function I {ml, m|) in 
(25), or Px/{mi + 1712) < 1 for the whole integral (26). Nevertheless, this is not an essential 
drawback. To get rid of it, one can apply the linear transformation (Ref. 1 p. 30, Eq. (83), 
Ref. 20 p. 218, Eq. (8.3.3)) 

/ 7/7/ — 7? \ 

Fi(a; b, b'; c; u, v) = {1 - uY'^Fi (a;c-b- b' , b'; c; -, ) (27) 

V u — 1 u — 1/ 

to the Appell function Fi in (25). This provides the necessary analytic continuation for 
this function, very similar to that given by the linear transformation (13) for the Gauss 

hypergeometric function: The first argument of Fi behaves in exactly the same way as the 
variable z in (13). Hence we obtain, for any mi, m2 > 0, 

i{mlml) = 2^ [l + (mi + m2)']~^ (28) 

TP f £ -i £^ 13 1 l + (mi-m2)^ 

X -f^i 7:; 1 



2' 2' 2 2'2'l + (mi+m2)2' l + (mi + m2)2^ 

The convergence range of this double series expansion now coincides with the initial validity 
domain of the original integral (1). 

Although with a reduced domain of convergence, the compact expressions (25) or (26) can 
be used in further calculations, where the needed analytic continuation is done at some later 
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stage or is provided by the extended convergence region of final results. Examples of this 
kind will be given in the following section. Also, the Appell function Fi is presently a built-in 
function in the Mathematical^, and its analytic continuation is automatically provided by 
this program. 



V. SPECIAL CASES 



Our general results (17) and (25)- (26) can be checked by speciahzing to several integer 
space dimensions D. Here we encounter mathematical simplifications allowing us to get the 
expressions for /(mf, m|) in terms of elementary functions. Some of them have already been 
calculated by other means before. 

At D = 1, the integrals (1) or (5) are elementary, and a short calculation using 
Mathematical"^ yields 

The same result immediately follows from (25): With £ = 3, the Appell function Fi simply 
decouples there into two geometric progressions. 

At -D = 2, the integral i{m\^m^) has been calculated in Ref. 9, in the context of the 
1/A^ expansion for uniaxial {m — 1) Lifshitz points in three-dimensional {d — 3) space. In 
this special case we have s — 2, and we obtain from the equation (17) 

/K.'«i)U, = lf.(i4;|.i4.^). (30) 

This simplified Appell function reduces through the relation (Ref. 19, p. 102) 

''^("'^'^ + "-^'^' {l-m-Yy{l-X){l-Y) ) 

= (1 - YY (a, 6; 1 + a - 6; -X^^) (31) 

to a Gaussian hypergeometric function with parameters 1,1/2; 3/2 and positive argument 
w — (mi — m2)^)/(l -I- (mi -\-m2Y) < 1. It is known (e.g. Ref. 18, p. 556) to represent 
a logarithmic function. Thus we obtain 



ft 2 2^l 2 l + y/^ 2 1 + m^ + mj + VA 

/(mi , mn) „ „ = —= In ^ = In . (32) 

V 1' 2)\D=2 ^ 1-V^ 2mim2 ^ ' 
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This special result comes as well from (25) through the reduction formula (Ref. 19, p. 79) 

Fi{a;b,b';b + b';u,v) = {1 - v)-%F,{a,b;b + b'; (33) 

followed by the linear transformation (13) for the resulting Gauss hypergeometric function. 
The latter transformation converts the argument z = {u ~ v)/{l — v) = —(1 + (mi — 
^2)^) / {4:17111712), which blows up for small values of nii and m2, again to the safe combination 
w < 1. This is a counterpart of the analytical continuation carried out before through the 
transformation (27) for the Appell function Fi. 

The last expression in (32) agrees with the integrand of 7(1, q) from equation (70) of Ref. 9 
after the shift of the integration variable via q' ^ q' — q/2 and identifications mi = {q' — q/2)'^ 
and 1712 = (q'' + 9/2)^. After some work, it can also be seen to be equivalent with the formula 
(4.3) of Ref. 11. 

At D — 3, the dimensional parameter s — 1, and in the Appell function Fi from (25) one 
of the nominator parameters vanishes. Thus, the function Fi reduces to an ordinary series 
expansion, and we get 

/(mi , ^71.2) |^^3 = 2 arctan(mi + m2)~^ (34) 

in agreement with Ref. 25, Eq. (A. 2)^6 and Ref. 11, Eq. (4.6). 

When D = 4, the value of e vanishes, and both the equations (17) and (25) yield the trivial 
result /(mi , ml)\j-)^^ — 1. Now, of interest is the first-order term of a small-£ expansion of 
/(mi, m|) taken eit D — 4 — e. Its calculation will yield the essential non-trivial contribution 
to the finite part of the whole integral I{ml,mf) given by equations (5)-(6). This term has 
been calculated previously by BDS^ using the Feynman parametrization. 

Expanding the result (17) in small e we obtain (see Appendix A) 

I{ml,ml) = l + £ + -VAln 1 1 ^ + - m^ - m^ In — | 

- |ln(mX) + 0(£') (35) 

in agreement with Refs. 2 and 11. 

The limiting case D — is somewhat special. It can be well illustrated by turning to 
the "double" integral (3) stemming from the Lifshitz-point theory. Let us accept again that 
D — d — m while the d-dimensional space is split into mutually complementing m- and D- 
dimensional subspaces. Then, by inspecting (3), it becomes evident that when the dimension 
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D shrinks to zero, we must remain with an m dimensional integral over q where no trace 
of D dimensional integration remains. This implies that in the limit — > 0, the result of 
the integration over p must yield its integrand aX, p = = Physically relevant limiting 
regimes of this kind have been considered in Refs. 8,9. For the integral I{px'iTn\,m'^) from 
(1) this means that 

I{Px]m\,ml)\^^^^-^. (36) 

This is easily reproduced from equation (26). At Px — 0, the first argument of Fi vanishes, 
and it reduces to a Gauss hypergeometric function yielding 

I{Px = 0;m?,m^) = (47r)-^+-/^r(i) / ~ 

\2J Z — e mi — 1712 

in agreement, up to the normalization of the integral, with Ref. 10. By setting e — 4 here, 
we obtain (36). Note that this correct hmit could not be obtained by setting D — directly 
in the scaling function (25). There, no vanishing of the first argument in Fi occurs, which 
was achieved by restoring the external momentum dependence in (26). 



VI. REDUCTION RELATIONS FOR THE APPELL FUNCTIONS 

By equating the results (17), (25), or (28) for the integral /(mf,m|) with its BDS ex- 
pression (7) we obtain apparently new reduction relations for the involved Appell functions. 
Thus, for the function F4 we get 

/ 1.3 1. 1 {n4^m^\ _ T^O:^^ .37^ 
^4"' 2' 2' 2' A' A r{2-2ay'^ ^^^^ 

+ 2(1^^0 j 



^K^J — w^) — ^^^^"'2-"' ) 

where a parameter a stands in place of 2 — D/2 — e/2. The relation (37) was derived for 

arbitrary real positive masses m\ and m\. This is in conformity with the standard definition 
of the convergence region of the Appell function F4, (1 + |m| — m^|)/\/A < 1, as discussed 



m 



2\2 



at the end of the section III. Once again we recall that A = 1 + 2{m\ + m^) + {ml 

The parameter a was originally constrained to the interval < a < 2 by physical appli- 
cability range of the integral (5). In (37), this limitation can be removed. The relation (37) 
holds for negative values of a, where the original integrals are divergent in the ultraviolet, or 
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for q; > 2 where the integrals would be negative dimensional. We have numerical evidence 
that (37) is valid for integer numbers a > 2 where the hypergeometric functions become 
singular and one might beware of possible exceptions. This is in agreement with the smooth 
dependence of F4 on a that indicates no peculiarities for positive integers. Moreover, in the 
relation (37), the parameter a can be considered to extend to complex values, along with 
the variables m\ and m|. 

Eliminating the parameters m\ and m| in favor of variables that appear in the function 
F4 we express the last relation in a standard fashion, 

-^(44.?-) --"g^;^ (38) 

2(1- a) \ v^^ + ^+l/ 

2(1- a) \ V^^-y^+l/ 
Comparing the equations (25) and (7) we write down a similar reduction relation for the 
Appell function Fi, 

F,(<.;l,a-i;|-/,<;^) =4-"^;|^3^(p=)-'+">sj-« (39) 

where 

s= {l + 'P^){'P^ + q^)- 

By writing in (39) the variables of Fi as —p^ and we stress that this function is even in 
p and q. So is the combination on the right, too. But here, each of terms with a Gauss 
hypergeometric function contains an explicit dependence on q and is not even in this variable. 
All the odd powers of q must cancel in the whole combination through the symmetrization 
/(g) + f{—q) present here. Its role is analogous to that of (2), implemented in (7) and (37). 

The Appell function Fi from the relation (39) exists for any real or complex values of 
the parameter a. It is absolutely convergent for p^ and < 1. These variables can also be 
both real and complex. The convergence region of the function Fi can be extended by using 
the transformation formulas like that in (27). A further reduction relation for the Appell 
function can be written in a similar way by using the formula (28). 
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VII. CONCLUDING REMARKS 



In this paper we have discussed the functional form of a standard but non-trivial one-loop 
Feynman integral (1) with non- vanishing masses (mf and m|) and external momentum in D 
spatial dimensions. Previously the results for this Feynman integral have been obtained by 
Berends, Davydychev, and Smirnov^ and reproduced later^^"^^ in several equivalent forms. 
The BDS result contained a linear combination of two Gauss hypergeometric functions 
neither of which was symmetric with respect to the interchange of the masses m\ and m\. 
The correct symmetric part of this combination had to remain after very complicated and 
obscured cancelations of non-symmetric terms. 

Here we represent the new explicit results for the Feynman integral (1) in terms of Appell 
functions, manifestly symmetric with respect to the masses m\ and m^. The advantage of 
our results over other representations is that they are expressed by a single function with 
completely symmetric and rather simple arguments. In equation (17) the result is given in 
terms of F4, "perhaps the most intractable Appell function" A further transformation of 
this function led us to elegant results (25)- (26) and (28) expressed in terms of the Appell 
function with extremely simple and symmetric arguments. 

The knowledge of new functional forms of the integral (1) gave us the possibility to derive 
the equalities relating the involved Appell functions to previously known combinations of 
Gauss hypergeometric functions. The way of producing the new mathematical identities 
by evaluating some integrals in different ways is not new. It was successfully employed by 
Inayat-Hussain^^'^° (see also Ref. 31) while considering certain Feynman integrals arising 
in phase transition theory, or by Srivastsva, Glasser, and Adamchik^^ while studying the 
Riemann Zeta function. 

Apart from equations (37)-(39) for F4 and Fi, similar reduction relations can be derived 
also for the Appell functions F3 and F2 since Fi can always be expressed in terms of one of 
these latter functions^ Of special interest would be the relations involving F2 in view of 
attention attracted by this function in the recent mathematical literature^^"^^. 
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APPENDIX A: THE EPSILON EXPANSION OF I {ml, ml) 

Let us consider the equation (17) at small e. For a while, we denote the arguments 1/A 
and 4ml/ A Appell function F4 as x and y, respectively. We use the definition (16) of 

F4 and split there the A; = term. Denoting it by ctq we have 

ao = E^^" = l + ^E- + ^(^') = l-|Ml-l/) + 0(s^) (Al) 

n>Q n>l 

since 

(a)„ = a(n-l)! + 0(a^) if n > 1. (A2) 

After a rearrangement of Pochhammer symbols via {a)n+k — (o)fe(a + k)n, the rest of the 
double series expansion can be written as 

s/2)k{l/2)k x" 



fe>l 



Again we use the property (A2) of {s/2)k to obtain 

= I E 7^ T V2 + k- 1/2; y) + 0{s'). 

k>l ^ ' 

Here, the hypergeometric function reduces to an algebraic function (Ref. 15, p. 34, Ref. 16, 
p. 461, Eq. 106) through the same relation that was used in (10). Thus we get 

4 1: wl ? [(^ + + (1 - vsn + o(e^). 



This formula implies that, if we introduce the series expansion 

^ (l/2)fc 
\ (3/2), k ' 



fe> 

the function a\ is given by 

= f K^+) + <^M+0(£^) where X^^ (A4) 
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By shifting the summation index in (A3), we can write a{z) as a a generahzed hypergeometric 
function 3F2: 

= ! 3^,(1 1,1; I 2;.). 
This can be expressed via (Ref. 16, p. 519, Eq. 366) in terms of logarithmic functions: 



a{z) = 2 - ^ In i±:^ _ ln(l - z) . 
y z 1 — \J z 

The same result could be reached by noticing that the derivative of o(z) is given by 



(A5) 



(1/2) 
(3/2) 



'2' 



' lnl±vS-l 



2^z 1-x/i 



where the last equality follows from Ref. 16, p. 477, Eq. 157, and integrating back o'{z) 
with respect to z. 

The summary of the above calculation is that the first-order e expansion of the Appell 
function F4 from (17) reads 



^/el31 \ e,, . e ( x \ / x \ 

2' 2' 2' ^' = ' - 2 - + 4 Httttt^J + Htw^J 



(1 + V^)^' 

where the function (t{z) is given explicitly by (A5). 

Finally, the small-£ expansion of the complete formula (17) appears in the equation (35) 
of the main text. 
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